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Abstract 



Chiral symmetry breaking due to instanton-produced fermion zero modes 
in the confining vacuum is considered. Zero modes provide 'tHooft-type de- 
terminant quark interaction, which is bosonized by introduction of auxiliary 
fields. For three flavours this procedure becomes nontrivial and a method is 
suggested which allows to derive effective Lagrangian for massless Nambu- 
Goldstone modes. 



1 



1 Introduction 



Recently effective action has been obtained for a theory with confining back- 
ground superimposed on instantons [1]. This is an extension of earlier investi- 
gations of the pure instanton vacuum [2-4]. The aim of [1] and of the present 
letter is to formulate a model of the QCD vacuum with properties of confine- 
ment and chiral symmetry breaking (CSB) at the same time. However in [1] 
only one- and two- flavour (Nf =1,2) cases have been considered. It is known 
that bosonization becomes really nontrivial for Nf > 2 [5]. In the present paper 
we perform explicit path integral bosonization of quark dynamics outlined in 
[1] for Nf = 3 case. 

First we briefly remind the basic features of the QCD model developed in [1]. 
The driving idea is to consider chiral and confinement effects as interconnected. 
To this end the QCD vacuum is assumed to consist of instantons plus confining 
configurations. Thus the starting point is the ansatz for the vacuum 

N 
i=l 

where A 1 ^ is the field of z-th (anti)instanton (N is the total number of instantons 
and antiinstantons in the 4-volume V) , B is the background field which ensures 
confinement (i.e. correlators F^(B) yield nonzero srtring tension). Instantons 
through fermion zero modes give rise to CSB. 

Confining background on one hand modifies the instanton density [6], and 
on the other it interplays with instantons and hence modifies chiral effects [1]. 

As shown in [1] the gauge invariant partition function of quarks in the field 
(1.1) in the limit of zero quark masses can be written in a form similar to that 
of the pure instanton gas [4] : 

Z = f Di/j Di/j + D fi(B) J dg+dg-expW , (1.2) 

W = -N + lng+ - N_lng- + f d 4 u{ipj(iD)ip f + g+detJ+ + g_detJ-} , (1.3) 

where the action of the gluonic field is included in D/i(B), the integration over 
g + and e/_ is a remnant of the inverse Laplace transform introduced in order to 
rewrite the partition function in the exponential form [4], / = 1, 2, Nf is the 
flavour index, N + and N_ were defined in Eq.(l.l), D = D(B^) is the covariant 
derivative corresponding to the field B^, the Nf x Nf matrices J± are related 
to 2Nf fermion vertices Y± via the relations 

Y± = J d A udetJ ± {u) , (1.4) 
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(J±(u))fg = J dxdyipj(x)-(l T ls)K(x, y, u)^ g {y) , (1.5) 

K(x, y, u) = iD(p{x — u)<$>(x, u, y)(p (y — u)iD , (1-6) 
where (p{x) is the instanton zero-mode solution: 

= ~ ( 2 P 2\3/2 "7=1 ( 1 - 7 ) 

7r (x 2 + p i yi i VX 2 

and $>(x,u,y) = $>(x,u)$(u,y) is a product of parallel transporters 

rx 

<S>(x,u) = Pexp (ig / B^dz^). (1.8) 

Note that J± and W are gauge invariant due to the presence of long derivatives 
D and parallel transporters in (1.6). 

Expressions (1.2-1.3) for the partition function are the starting point for 
the rest of the paper. Our goal is to derive an effective action by integrating 
over the fermionic degrees of freedom. The most efficient way to do this is 
to introduce auxiliary fields (bosonization), then perform the integration over 
fermions exactly, and then apply the stationary-phase approximation to the 
integral over the auxiliary field. As a result one obtains the effective action for 
Nambu-Goldstone modes. 



2 The case of two flavours 

The usual path integral bosonization is based on the Hubbard-Stratonovich 
transformation [7] . For the quark dynamics under consideration this procedure 
has been performed in [4] and [1]. However it faces problems for Nf > 2, and 
more sophisticated manipulations are needed. In this section we shall reconsider 
the Nf = 2 case (see Ref.[l]) within the framework of the approach applicable 
for arbitrary Nf. The basics of the bosonization techniques which will be used 
here have been developed in [8,9]. 

Consider the partition function (1.2-1.3) for Nf = 2. The general relation, 
which holds for any 2x2 matrix, allows to write 

det J± = \[{trJ ± f - trJ 2 ± ] . (2.1) 

Expanding J± over a complete set of three Pauli matrices r a (a = 1,2,3) plus a 
matrix r° = —il yields 

J±=tr a c^ C J =1)2j3 = \tr{r a J ± ) , cf = ~tr(r°J ± ). (2.2) 
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Then by virtue of (2.1) we get 
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det J ± = - Y. (ci? ■ (2.3) 



a=0 

Next we split out from the partition function Z given by Eqs. (1.2-1.3) the 
following part: 

Z = J D^D^ + exp{ J d 4 u[ipj{iD)ip f + g+detJ+{u) + g^detJ^u)}} (2.4) 

= J DipDip + expW. 

This expression can be written in the following equivalent form (index a runs 
from to 3): 

Z = j DifjDi/j + n Dai II Da-5(ai - g + c+)5{cj- - g^-)exp(W) = 

J a a 

= / Dt/jD^ n Da+ a n Do' R DM+ ft DM a exp(W) x 

J a a a a 

x exp{i J d A u[M+(o+ - g + c + a ) + M a (o~ - g. C -)]} . (2.5) 

At the next step we integrate out in (2.5) the quark fields tp and i.e. 
perform the integration 

Z^ = J D^ + exp{i j d 4 u[^fDiJf - g + M+c+ - g-M^c')}} , (2.6) 

where according to (2.2) and (1.5) we have 

1 a fg 

= J d 4 xd 4 y^(x)[^M ± (u)( ) ^)K(x,y,u)}tf(y) , (2.7) 
with M ± (u) ee {M± >0 (u)r a>0 , - M ± r }. The integration yields: 

Z^ = -exp(ln detX) (2.8) 

/ -i(D A - i&D) ^g+M+K \ 
X= ' . (2.9) 

v \g_M~K -i{D 4 + iBD) , 

Next one performs integration over the fields 

Z a = f [] Do+ J] Do~ exp{ f d A u[g+detJ+ + g-detJ- + 1M+0+ + zM~ct~]} . 

(2.10) 



Making use of (2.3) and (2.5) one has 

g± detJ ± = - { -^- , (2.11) 
9± 

where summation over a is implied. Therefore integral (2.10) is Gaussian with 
the result 

Za = exp{-± J d\[g + (M+f + 9 -{M-f]} . (2.12) 
With the identification of the fields with L a and R a of Refs. [1,4] 

La = ^M+, Ra = \M~ . (2.13) 

the final results for fully bosonic effective action reads 

exp{— S(L, R)} = J dg + dg-exp{— N + lng + — N-lng-+ 

+ J d A u[g+L 2 a {u) + g-R 2 a {u)] + In det X , (2.14) 

/ -i(D A - i&D) g+LK \ 
X= . (2.15) 

v g RK ~i(D A + iaD) j 

For further purposes it is instructive to rewrite the term In det X in (2.14- 
2.15) in a slightly different form. In (2.14-2.15) the determinant in Dirac space 
is written down explicitly. Instead let this step be implied in det X notation. 
Then making use of the identity In det X = tr In X one can recast In det X in 
(2.14) into the form 

IndetX = tr\n{ib + g + {^—^)LK + g-^-^)RK} , (2.16) 

This form will be used in Section 4 in order to define the chiral effective 
Lagrangian. 

We have retrieved the results obtained in [1]. The reader is referred to [1] 
for a discussion of the physics behind Eqs. (2.14-2.16). For Nf = 2 our present 
derivation is more cumbersome than the Hubbard-Stratonovich transformation 
used in [1]. However our goal was to display the algorithm applicable for Nf > 2 
case. 
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3 Bosonization of Three Flavours Dynamics 



Again our starting point is the partition function given by Eqs.(1.2-1.3), and 
our task is the integration over fermionic degrees of freedom. 
For the case Nf = 3 relations (2.1-2.3) are replaced by 

det J± = \{trJ ± f - htrJ ± )(trJl) + \trJl , (3.1) 

D 2 6 

J±=J2 A a cJ , c a>0 = \tr\ a J ± , co = \tr{\ Q J±) , (3.2) 

a=0 ^ o 

detJ ± = (4f~yj:(c±) 2 + \ E tr(\ a \ b \ d )c±c±4 ■ (3-3) 

z a=l <J a,6,d=l 

Here A a , a = 1, 2, .., 8 are the SU(3) Gell-Mann matrices, A a=0 = 1. 

Next one again proceeds as in equations (2.4-2.6) with the only difference 
that the index a now runs from to 8. In (2.7) the r a matrices are replaced by 
A a , so that now M ± {u) = M±{u)X a . In line with (2.8-2.9) and (2.16) integration 
over Dil)Dil) + yields 

Z^ = exp{trHib + g + ^^)LK + g^^^)RK]} . (3.4) 

Here L = \ M+ = L a X a , R = \M~ = ,R a A a , so that (3.4) contains 18 bosonic 
fields L a , R a , a = 0, 1, 8. 

The next step is the integration over and that is where the real difference 
from the Nf = 2 case starts. The integral analogous to (2.10) reads 

1 f\ 8 1 8 

Za = I II^aexp{- / d A u[--{al - -a E °\ + ^ E tr(X a X b X d )a a a b a d ) -(3.5) 

a J g Z a= i 6 fl; 6 j( f=l 

-* E M a a a }} = / l[Da a exp{S(a a )}. 

a=0 J a 

In (3.5) we have suppressed the (±) indices, so that additional summation over 
(±) is implied. The Nf = 2 integral (2.10) was Gaussian and therefore it 
could be calculated both exactly and by stationary phase method yielding the 
same results. Exact analytical calculation of (3.5) is not possible. Therefore one 
has to resort to the stationary phase, or the steepest descent methods. The 
corresponding conditions which define the a a = f({M a }) solutions have the 
form 

<)S ~ 3^,i ^ E ^ • '/r-U) - o 



5a 



0=1 
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5S 1 8 

-— = -5a a a + - E tr{X a X b X d )a b a d + zp 2 M a = , a = 1, 2, 8 . (3.6) 

The principal question is the existence of a contour along which the integral 
(3.5) converges, so that the steepest descent method yields reliable results. As 
a toy model consider the case <Jo 0, <j a = 0, a = 1, 2, 8 (for Nf = 2 similar 
situation was investigated in [1,4]). Then (3.5) reduces to 

3 

Z a = J Daexp{- J d 4 u(-^ - iMa)} (3.7) 
with a — do- Replacing M = —2iL, g = in, we get 

Z a = J Daexp{ J d A u(2La - ^)} . (3.8) 

This is an Airy type integral [10]. In the Appendix we show how to choose 
the correct contour and how to calculate this integral by the steepest descent 
method. Up to unimportant overall factor the result is 



k J d 4 uL 3 / 2 (u)} . (3.9) 



This result is reminisced of the stationary phase condition for Landau magneti- 
sation function of the infinite range Ising model [7]. 

4 The gap equation 

For the toy model under consideration the partition function (1.2-1.3) reads 



Zoc J DL J dKexp{-N In k-^^k J d 4 uL 3/2 + trln(-D 2 + k 2 L 2 K 2 )} . (4.1) 
The steepest descent integration yields 

- 2 #^ 1/2 + 2k2l tr U^W = • 



2 



N 4 



2_o/o _ ,o . / K 



7 -^V^ + ^ __ ) = o, (4.2) 
From (4.2) we find 



2r2, K l 3 



With the identification for the chiral mass operator m = kLK, and performing 
on the l.h.s. of (4.3) trivial summation over flavour indices one gets the same 
form of the gap equation as for Nf = 1, 2 in [1] 

m 2 , N /A A ^ 

= y • (4-4) 

5 The effective chiral lagrangian 

We turn now to the final topic, namely to the constuction of the effective 
chiral Lagrangian for Nf = 3. Such a Lagrangian corresponds to an octet 
of Goldstone bosons 7T°, n ± , K ± , K°, K°, j] which are coupled by the fermionic 
loop. The starting point is Eq.(3.4) in which we parametrize bosonic fields L 
and R in the following form [4] 

L = L a X a = ip( 1 + ^ + ^ )UV, R = R a X a = ip( 1 + ( 2 + fl )VU + , a = 0, 1, 8, 

1 1 (5 ' 1} 

U = exp{^7r fe A fe } , V = exp{^ie b \ b } , 6=1,2, 8. (5.2) 

This is a general parametrization which expresses 18 bosonic fields L a , R a in 
terms of 18 fields £, fi, 7Tb, £&(note that in (5.1) summation is from to 8, while 
in (5.2) from 1 to 8). Insertion of (5.1)-(5.2) into (3.4) yields 

= exp{tr \n[W + z(l + £ + /i)UVE+ + i(l + £ - fi)VU + E-}} , (5.3) 

where E ± = EQ^-) , E = \gpK, and the operator K is defined by Eq.(1.6). 
In order to get the effective chiral Lagrangian one has to integrate (5.3) over 
all meson fields exept for the Goldstone octet Using the steepest descent in 
the vicinity of the vacuum sadde-point [1,4] we get 

Seff(n) = - In / D^D/iDeZ^ = 

= tr\n{iD + iEU h ) , (5.4) 

where 

U 5 = t/(^p) + CZ+f^p) = expi^A-Ts} . (5.5) 

The form of the effective action (5.4) coincides with that found in [1] for Nf = 2. 
In absence of background gluonic field the form (5.4) reduces to the effective 
action found in [4] for the pure instanton case (see also discussion in [11]). 
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V.A.Novikov for helpful discussion and to D.Kuzmenko for his contribution 
to the Appendix. 

The final note is that the authors became aware of Ref. [5] after completing 
the present paper. Although our approach differs from that of [5], the results 
we got for Nf = 3 fit into the sheme outlined in [5] for arbitrary Nf. 

APPENDIX 

In the Appendix we show how the result (3.9) emerges. Consider the integral 
(3.8): 

3 

Z = J Daexp{ J d 4 u(2La - ?-)} . (A.l) 
We introduce a new variable w and a new parameter t according to 



w = 



'2k 2 L\ _1/2 (2\ 1 ' 2 V2 

a, t = 2(-\ kL I . (A.2) 



Then (A.l) takes the form (up to a Jacobian which is not essential for our 
present purposes) 



Z = j Dwexp{ j d 4 ut{w - — )} , (A.3) 



3 

which is the standard representation of the Airy function [10]. The leading 
contribution (3.9) is then obtained by the steepest descent method. The integral 
has two stationary points w = ±1, from which only w = — 1 contributes. The 
integration contour can be choosen along the line Rew = — 1. 
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